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A.  DEVELOPMENT  OF  A  LIMITED-AREA  NUMERICAL  WEATHER  PREDIC¬ 
TION  MODEL 

I.  INTRODUCTION 

Numerical  Weather  Prediction  (NWP)  models  of  various 
forms  have  been  used  operationally  in  the  United  States  for 
about  25  years.  Their  increased  sophistication  and  complex¬ 
ity  have  paralleled  that  of  the  high-speed  computers  on 
which  they  are  implemented.  They  have  developed  from  the 
simple  equivalent  barotropic  model,  which  used  various 
simplifying  assumptions  to  predict  the  height  of  the  500  mb 
pressure  surface,  to  the  seven-layer  primitive  equation 
model  currently  in  use  by  the  National  Meteorological  Cen¬ 
ter  ( NMC ) . 

Models  have  developed  in  complexity  as  a  result  of 
attempts  to  account  more  adequately  for  the  various  physical 
processes  which  signif icantly  affect  weather  changes.  One 
way  is  by  numerical  integration  of  the  so-called  "primitive 
equations";  that  is,  the  spatial  and  temporal  mathematical 
relationships  in  their  basic  form  with  as  few  simplifying 
assumptions  as  possible,  which  describe  atmospheric  behavior. 
Because  of  computational  storage  and  time  requirements,  the 
equations  are  integrated  in  time  over  a  spatial  network  of 
grid  points  of  limited  resolution  in  both  the  horizontal 
and  vertical.  Since,  however,  many  important  physical  pro¬ 
cesses  have  spatial  scales  smaller  than  the  computational 
grid  interval,  they  cannot  be  properly  accounted  for  by  the 
system  of  approximate  equations.  In  addition,  the  larger 
the  grid  spacing,  the  greater  the  error  associated  with  the 
finite  difference  techniques  used  in  solving  the  nonlinear 
partial  differential  primitive  equations. 

On  the  other  hand,  the  observed  data  used  to  provide 
initial  and  boundary  values  for  the  model  are  available  only 
at  a  density  prescribed  by  the  density  of  the  observati onal 
network.  Some  information  may  be  added  or  lost  by  interpo¬ 
lation  to  accommodate  a  grid  of  finer  resolution.  There- 


of  the  hydrostatic  meteorological  equations,  given  here  in 
Cartesian  form  on  a  conformal  projection  of  the  earth: 


For  this  system  of  equations  to  comprise  a  closed 
set,  the  following  accompanying  diagnostic  equations  were 
included  in  the  model: 


In  these  equations  x  and  y  are  the  two  horizontal  Cartes iar 
coordinator  on  the  projection;  t  time,  u  and  v  the  x  and  y 
component  s  of  velocity  true  on  the  earth;  6  potential  tem¬ 
po  ra  f urn ;  ♦  -  gz ,  geopotent ial  above  the  projection  where 
i'  is  a  ■■■  lerat  ion  due  to  gravity  (taker,  as  constant)  and 
■  :rht  above  sea  level;  p  pressure;  1,  1,000  mb; 

*?.■  nr  ■  ■  ■  t:s  ■  ant.  for-  dry  air;  c  the  specific  heat  of  air 
a  *  ■"*  a:.4  score;  m  =  (1  +  sin  (•0°)/{  1  +  sin  $),  where 

0  ■  ■  *  f  :  :  *  :  ‘’iif  -,  th>  rr.ap  factor-  of  the  polar'  stereo- 


graphic  projection  with  60°N  as  the  standard  latitude;  f  the 
Coriolis  parameter;  and  W  the  precipitable  water.  Fx  and  F 
are  the  x  and  y  components  of  acceleration  due  to  friction 
and  turbulence;  H  represents  diabatic  heat  influences;  and 
C  the  rate  of  gain  or  loss  of  water  vapor  due  to  evaporation 
or  condensation.  Vertical  velocity  is  o  =  do/ dt. 

Vertical  coordinate  a  is  expressed  in  the  general  form 


P  ~ 

Pl  -  Pu 


(9) 


where  is  the  pressure  on  a  given  quasi-horizontal  surface 
bounding  an  atmospheric  layer  above,  p^  the  pressure  on  the 
surface  bounding  the  layer  below,  and  p  the  pressure  at  a 
point  in  question  within  that  layer.  Coordinate  a  is  de¬ 
fined  differently  in  each  of  the  four  domains,  as  illustrated 
in  Fig.  A-l ,  which  depicts  the  vertical  structure  of  the 
model . 

The  area  chosen  for  consideration  in  this  study  is  a 
segment  of  the  NMC  octagonal  grid  which  completely  includes 
the  contiguous  United  States.  The  coarsest  resolution  used 
in  this  study  is  equal  to  twice  that  of  the  NMC  octagonal 
grid,  with  grid  points  190. 5  km  apart  at  60°N.  A  grid  of 
33  x  33  grid  points  was  used  over  which  computations  were 
carried  out,  with  the  outermost  row  designated  as  the  boun¬ 
dary.  The  model  was  initially  designed  to  consider  only 
this  resolution,  but  is  being  modified  +o  accommodate  double 
and  four  times  this  resolution  in  order  to  accomplish  the 
primary  goal  of  the  study. 

Numerical  solution  of  the  model  equations  begins  with 
the  specification  of  initial  values  of  prognostic  variables 
u,  v,  0  ,  3p/3o  ,  and  W.  Values  of  variables  u,  v,  and  0 
must  be  specified  initially  for  all  seven  levels,  but  0  is 
held  constant  in  time  and  space  in  the  uppermost  or  cap 
layer.  Holding  8  constant  in  this  way  facilitates  use  of 
the  nap  layer  for  strictly  computational  purposes  -  the 
valuer  of  u,  v,  and  ap/so  forecast  in  this  domain  are  used 
only  to  obtain  values  for  the  variables  in  the  other  domains. 
Initial  values  for  W  must  be  available  at  all  grid  points 
in  the  lowest  three  layers  only,  finer  it  is  taken  to  be 
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.ure  A-l .  Schematic  diagram  of  LFM  vertical  struc 


9 


zero  everywhere  above.  W  is  obtained  in  each  -layer  from 
the  initially  specified  specific  humidity  q  by  means  of  the 
relat i onsh i p 

r°  L 

W  =  JJq  do  (10) 


where 
quest i on 


and  o ^  represent  top  and  bottom  of  the  layer 
Fip.  A-P  illustrates  the*  relative  locations 


i  n 
of 


the  prognostic  and  diagnostic  variable  values  in  each  layer. 

Once  these  values  arc  obtained  the  model  determines 
the  pressures  p  at  all  prid  points  on  the  interfaces  between 
layers  by  vertically  intopratinp ,  usinp  values  of  3  p/ 3o  • 
S.ince3p/3o  is  independent  of  hoi, "'lit  in  each  domain,  pres¬ 
sure  at  each  o  interface  in  each  domain  is.  fiver,  by  the 


express i on 


p(  o) 


p  ( o  =  0  ) 


til) 


With  pressure  at  the  top  of  the  model  atmosphen  set.  .  qm 
to  zero,  this  expression  is  used  to  evaluate  the  \  r<  sou:- 
all  prid  points  (with  pg  =  p1  -  p„  -  ‘>0  mb  ev.  rywh<  ••• 
boundary  layer).  From  these  values  the  model  o!  s  . 
for  w  from  (?),  so  that  (e)  may  be  integrated,  num.  s'  .  . 
upward  to  obtain  the  values  of  $  at  all  layer  :r.:<  s' 

At  this  point,  (8)  is  evaluated  numerically  at  a..  s 


interfaces  except  for  the  materia  1  surfaces  desipna  t  <  d  by 
k  =  1  ,  8,  7,  and  8  (see  Kip.  A- 1  )  where  i  t  is  assur.t  d  *  hat 
a-  0.  These  assumptions  also  serve  as  boundary  .orditlo: 
in  so.lv  inf  (8).  In  addition  to  these  boundary  conditions, 
tho  boundary  condition  for b  ^  at  the  top  of  the  boundary 


layer 


rp  \  v  rp 


(l.  ) 


where  Up  and  Vp  are  the  boundary  layer  velocity  components, 
is  obtained  by  i n t epra t i np  tho  boundary  layer  version  of 
the  eonfir.nify  equation  (3)  throupli  the  depth  of  the  l  curd 
■my  layer  and  by  observinp  the  requirement  of  mass  balance 
i  ::  tic  t  ( ■  1  a  1  *  tu  i  pc. '  ph  e  r ,  ■  . 

Te  el  !  a  i  ! .  'hi  p  r  ■  ■ .  j  i  e  !  ,  ■  d  Villi.  id'  W  foe  M,,.  love.'!  *  (if. 


i  I 


1  0 


tssuttied  t  ha  ! 


lea  i 


each  grid  point  is  a  linear  function  of  pressure  through 
all  three  layers: 

«  g  a  W  (13) 

where  a  is  the  function  of  0  that  will  determine  the  distri¬ 
bution  of  q.  When  ( 1 3 )  is  substituted  into  the  expression 
obtained  from  combining  the  conservation  of  specific  hu¬ 
midity  with  (3) ,  the  prognostic  equation  (5)  for  W  is  ob¬ 
tained.  Values  of  a(x,y,o,t)  in  each  layer  at  each  time 
step  are  obtained  by  means  of  the  expressions 

°B  '  aB  +  V B 

in  the  boundary  layer  and  two  tropospheric  layers  respec¬ 
tively,  where  a  and  b  are  arbitrary  functions  of  x,  y,  and 
t  only.  These  are  derived  from  four  simultaneous  equations 
(see  Shuman  and  Hovermale,  1968)  and  are  given  by 
brp  =  18/(2  +  3c)2,  bB  =  c2bT,  aT  =  -b^/3.  and  aB  =  2cbT/3 
where  c  =  (3p/3o)B/(  3p/)a) rj, . 

The  model  was  formulated  so  as  to  neglect  frictional 
influences,  diabatic  heating  effects,  and  changes  in  total 
water  vapor;  thus  Fx=F^=H=C=0.  This  was  done  so 
that  the  relationship  between  the  basic  dynamics  of  the 
model  and  model  resolution  might  be  investigated  at  first 
without  the  added  complication  of  these  more  complex  effects. 
In  addition,  water  vapor  concentration  has  been  set  ini¬ 
tially  to  zero  everywhere,  to  consider  the  simplified  case 
of  a  dry  atmosphere.  Values  of  Coriolis  parameter  and 
map  factor  were  calculated  for  all  grid  points  and  stored 
in  a  form  readily  accessible  by  the  model.  Eventually 
terrain  height,  sea  surface  temperatures,  and  drag  coeffi¬ 
cient  values  for  all  relevant  grid  points  will  be  made 
available  to  the  model.  As  the  model  is  now  formulated, 
the  only  one  of  these  parameters  necessary  is  terrain  height, 
which  is  set  equal  to  zero  at  all  ground  surface  grid 
points.  Thus,  at  this  point  the  model  has  determined  all 
values  needed  to  evaluate  the  right  sides  of  equations  (1) 

12 


through  (5)  as  they  are  currently  formulated. 

Terms  on  the  right  side  of  the  prognostic  equations 
were  coded  in  the  model  according  to  a  finite  difference 
scheme  devised  by  Shuman  and  Hovermale  (1968).  Since  dif¬ 
ferent  (albeit  perfectly  acceptable)  forms  of  the  finite 
difference  representations  of  the  terms  would  in  general 
yield  substantially  different  model  results,  great  care 
was  taken  to  calculate  the  terms  in  the  model  according  to 
the  methodology  of  Shuman  and  Hovermale.  When  the  finite 
difference  forms  of  the  right  side  terms  are  evaluated  at 
their  respective  grid  points  at  each  time  step,  the  model 
either  (1)  multiplies  their  sum  by  one  time  interval  in 
seconds  and  adds  this  to  the  initial  value  of  the  prognostic 
variable  in  question  for  the  first  time  step,  or  (2)  mul¬ 
tiplies  their  sum  by  twice  the  time  interval  and  adds  to 
the  value  of  the  prognostic  variable  at  the  previous  time 
step  for  successive  time  steps.  This  requires  storage  of 
two  arrays  of  values  of  each  of  the  prognostic  variables  at 
each  time  step,  in  addition  to  storage  of  the  values  of  the 
diagnostic  variables  at  the  present  time  step.  This  large 
data  storage  requirement  prohibited  simultaneous  storage  of 
all  values  necessary  for  calculation  of  all  variables  in 
core  memory  at  the  same  time.  Using  mass  storage  and  com¬ 
mon  array  storage  within  the  model  and  making  calculations 
for  each  prognostic  equation  in  separate  sections  of  the 
program  greatly  reduced  core  storage  requirements. 

Til.  ONE-LAYER  VERSION  OF  LFM 

To  elucidate  the  relationship  between  spatial  resolu¬ 
tion  and  accuracy  of  model  results,  it  is  desirable  to  inves 
tigate  progressively  more  complex  forms  of  the  model,  be¬ 
ginning  with  its  simplest  form.  It  was  believed  that  under¬ 
standing  the  effect  of  change  of  resolution  and  finite  dif¬ 
ference  forms  of  the  derivatives  at  simpler  levels  of  model 
complexity  would  facilitate  understanding  model  behavior  in 
its  more  complex  forms.  For  this  reason  the  three-dimension 
al  model  described  abovp  was  modified  to  include  logical 
switches  allowing  simulation  of  the  prognostic  variables 

1 1 


u,  v,  Q  ,  and  ap/Jto  in  one  atmospheric  layer  only.  This  is 
equivalent  in  theory  to  holding  vertical  velocity  equal  to 
zero  everywhere.  The  one-layer  model  atmosphere  is  con¬ 
fined  between  two  material  surfaces,  so  that  ♦  j  ,  ir^,  p^ 
and  ir0 ,  represent  the  values  of  these  variables  at 

the  top  and  bottom  a  surfaces  respectively,  with  =  1, 
o0  =  0 .  Prognostic  variables  u,  v,  and  9  are  considered  in¬ 
dependent  of  height  in  the  layer,  along  with  3 p/  3o  , 
which  is  also  height  independent  in  each  domain  in  the 
multi-layer  model. 

First  model  runs  were  thus  based  on  the  assumption  of 
a  dry,  adiabatic,  inviscid,  one-layer  atmosphere.  In  addi¬ 
tion  it  was  assumed  that  j  =  0  everywhere  (smooth  ground 
surface)  and  that  p0  =  0  everywhere  (top  of  the  layer  is 
top  of  the  atmosphere).  With  these  assumptions  the  prog¬ 
nostic  equations  (1)  -  (4)  reduce  to  the  form 


with  the  diagnostic  equations  =  c^eir,  a  (p)  -  p/p^  , 

3  p/3o  =  pj  ,  and  =  (pj/P)^cp,  The  modified  one-layer 
version  of  the  model  numerically  integrates  the  finite  dif¬ 
ference  form  of  these  equations  ever  the  33  x  33  grid  point 
grid  lattice  using  six-minute  time  steps.  Initial  and 
boundary  values  for  the  four  prognostic  variables  used  to 
obtain  the  numerical  solutions  to  these  equations  will  be 
discussed  below. 

Tn  addition  to  this  general  one- layer  version  of  the 
model,  two  simpler  versions  were  formulated  by  the  succes- 
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sive  incorporation  of  these  assumptions! 

(1)  Incompressible  motion  (assume  3p/do  =  1,000  mb 
in  time  and  space).  Model  equations  take  the  form 


where  ♦  n  =  c^O  . 


(2)  Incompressible,  inertial  motion  (assume  3p/3o 
=  1,000  mb  and  9  =  ?50°K  in  time  and  space).  In  this  case 
the  model  integrates  the  equations 


These  assumptions  were  individually  incorporated  into  the 
model  by  means  of  loci  cal  switches-  which  excluded  terms 
that  were  affected  by  the  respective  assumpt i ons . 

Initial  and  boundary  values  used  for  all  forms  of  the 
one-layer  version  of  I.I’M  were  derived  from  a  spatial , 
temporal  stream  function  which  is  a  solution  to  the  equa¬ 
tion  expressing  the  conservation  of  absolute  vortioity. 

This  equation  is  obtained  from  the  system  of  equations  des¬ 
cribing  two-dimensional,  i neompressi M e  motion  on  the  sur¬ 
face  of  a  spheres 
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(27) 


where  X  represents  longitude  and  <j>  latitude.  The  motion 
governed  by  these  equations  conserves  absolute  vorticity, 

K  +  f ,  where  £  =  (a  cos  0  )  ^  3  v/3X  -  3(u  cos  0  ) /34>  is  the 

relative  vorticity.  This  can  be  shown  by  cross  differentiat 
ing  (25)  and  (26)  and  subtracting  the  resulting  form  of  (25) 
from  that  of  (26)  and  by  observing  that  the  motion  field  is 
non-divergent  as  stated  by  (27).  The  result  is  the  expres¬ 
sion  for  conservation  of  absolute  vorticity: 


—  (  r  +  f)  =  -  - - — 

3t  ^  a  cos  0 


3X  ^  +  a  30  ^ 


+  f). 
(28) 


Because  the  motion  field  is  non-divergent  as  specified  by 

(27)  the  motion  field  V  can  be  expressed  in  the  form  of  a 

stream  function  0  such  that  V  =  k  x  V  0  ,  where  k  is  the 

s 

unit  vector  normal  to  the  spherical  surface  and 
V  o  =  (a  cos  0  )  i  3/3X  +  a  j  3/30  ,  with  i  and  j  repre¬ 
senting  unit  vectors  directed  along  latitude  and  longitude 

lines  respectively.  Thus  u  =  -a  *  30  /  30  ,  v  =  (a  cos  0  ) 

2 

30  /  3X  ,  and  £  =  V  o  0  .  When  these  expressions  are 

substituted  in  (28)  it  can  be  verified  that  a  stream  func¬ 
tion  of  the  general  form 


0  (0 


t)  =  A(t)F™  (  sin 


) 


im  X 


(29) 


is  a  solution  to  the  resulting  equation,  where  A  =  A^ 

exp  (  im2Qt/ n(n  +  1)}  in  which  A  is  an  arbitrary  initial 

u  - 1 

amplitude  and  Q  is  the  earth's  rotation  rate  in  sec  ,  and 
P™  (sin  0)  is  the  associated  Legendre  function  of  degree  n 
and  order  m.  This  travelling  wave  solution  to  the  vorticity 
equation  (28)  is  known  as  the  Rossby-Haurwi tz  Wave,  where  m 
and  n  specify  zonal  and  meridional  numbers  respectively. 
Initial  and  boundary  values  for  u  and  v  in  all  versions  of 
the  one-layer  model  were  derived  from  this  stream  function 
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for  the  case  m  =  2,  n  =  3,  and  were  transformed  from 
spherical  to  Cartesian  coordinates. 

In  all  but  the  case  of  incompressible,  inertial  mo¬ 
tion,  initial  and  boundary  values  of  potential  temperature  0 
were  required.  These  were  obtained  in  the  following  manner. 
It  was  first  determined  that  if  the  motion  dictated  by 
equations  (25)  -  (27)  was  to  be  confined  to  two  dimen¬ 
sions  (i.e.,  to  the  surface  of  the  sphere),  the  motion  field 
must  not  only  be  initially  non-divergent ,  but  that  the  di¬ 
vergence  tendency  must  also  be  zero  everywhere  in  the  motion 
field  in  order  to  conserve  mass.  This  additional  constraint 
is  used  with  equations  (25)  -  (27)  to  obtain  the  balance 
equation 

=  21 


V  2  * 

s 


[( 


a  cos  <}> 


[l  3c 

[a  3<J> 


if  J 1 

3$ 


anA/l  a 

3X  /\a 

2 


\U 


a  cos 


3  v 

31 A 


du 
3  <t> 


)] 


u  u  tan  6 
a 


+  1  3v  v  tan  ft 
a  34>  a 


(30) 


which  expresses  the  relationship  that  must  exist  between 
the  geopotential  field  and  the  motion  field  to  insure  that 
the  motion  field  is  instantaneously  and  temporally  non- 
divergent  (i.e.,  in  "balance"  with  the  geopotential  field 
at  all  times).  By  substituting  (29)  in  this  expression 
one  obtains  the  relationship  between  geopotential  field  and 
stream  function  with  which  it  is  in  balance.  Geopotential 
field  *  (  ♦  , X  ,  t)  that  satisfies  this  equation  was  expressed 
as  a  linear  combination  of  spherical  harmonics.  By  noticing 
from  (16)  that  non-divergent  motion  is  equivalent  to  holding 
3 p/3o  (and  thus  n^)  constant  everywhere  in  the  layer, 
potential  temperature  was  obtained  from  the  balanced  geo¬ 
potential  by  the  expression  0  =  $  2/cp,  where  *2 
=  2  ♦(♦  ,  A  ,  t)  from  the  assumption  that  ♦  in  the  middle  of 
the  layer  is  a  simple  average  of  the  values  of  geopotential 
at  the  top  and  bottom  of  the  layer.  In  this  way  initial 
and  boundary  values  of  9  (which  are  in  balance  with  initial 
and  boundary  values  of  u  and  v)  were  obtained.  Since  initial 
and  boundary  motions  are  non-divergent,  initial  and  boundary 
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values  for  9  p/9cr  were  held  constant  at  a  value  of  1,000  mb. 

While  balance  is  maintained  in  the  initiaj  and  boundary 
values  of  u,  v,  and  0  throughout  time,  no  such  guarantee  of 
instantaneous  and  temporal  non-divergence  exists  for  any  of 
the  forms  of  the  one-layer  model  equations.  The  balance 
equation  is  not  used  in  the  model  to  insure  such  a  balance. 

In  the  simplest  case  of  incompressible,  inertial  motion,  a 
divergent  residual  remains  when  the  right  side  of  (02)  and 
(21)  are  evaluated  using  the  non-divergent  initial  motion 
field.  In  the  other  two  cases  of  the  one-layer  model,  poten¬ 
tial  temperature  is  calculated  according  to  the  conservation 
of  potential  temperature  and  not  according  to  the  balance 
equation,  with  no  assurance  of  a  balance  between  <!>  and  u,v. 
Hence  results  from  any  form  of  the  one-layer  model  would 
not  in  general  be  non-divergent ,  while  the  initial  and 
boundary  motions  used  in  the  model  are  strictly  non-divergent 
at  all  times.  As  the  magnitude  of  the  divergent  component 
in  the  interior  motion  field  grows,  it  will  become  more  and 
more  unlike  the  non-divergent  boundary  motions,  and  disconti¬ 
nuities  may  arise  at  the  interface  between  the  interior  and 
boundaries  of  the  motion  field.  However,  the  balanced  mo¬ 
tion  and  temperature  fields  will  be  used  as  initial  and 
boundary  values  until  values  more  compatible  with  the  in¬ 
terior  can  be  found. 

A  smoothing  operator  was  also  incorporated  into  a 
separate  version  of  the  model.  It  is  applied  to  all  prog¬ 
nostic  variables  at  each  time  step,  and  is  identical  to 
that  used  by  NMC  as  described  by  Gerrity  and  Newell  (1976). 
The  smoother  is  a  2 5- point  operator  at  all  points  two  or 
more  grid  intervals  from  the  boundary  and  a  nine-point 
operator  at  points  only  one  interval  from  the  boundary. 

TV.  CURRENT  STATUS  OK  I  KM  M0DK1  DEV El  0  PM ENT  AND  I  IANS 

Some  preliminary  model  runs  have  been  conducted  with 
ail  three  versions  of  the  one-layer  model.  The  earliest 
runs  consisted  primarily  of  efforts  to  debug  the  model, 
with  only  the  most  recent  runs  yielding  results  which 
accurately  represent  the  dynamics  of  t  he  mode  1  equations. 

1  K 


Computer  rodps  are  currently  bo  i  ny  written  to  present 
results  in  yraphi  eal  form  so  t  hey  may  be  mo»-e  easily  in¬ 
terpreted  . 

Outputs  from  the  three  vnrsi  one  of  the  one-  layer  model 
will  be  analysed  and  compared  ayainst  each  other  to  deter¬ 
mine  the  effort,  o*'  each  sot  of  aesunv  1  ;  oils .  The  truncation 
error  of  its  numerical  finite  dl  f former  scheme  will  be 
analysed  by  numerieally  eompu  t  i  t:y  the  terms  on  the  right 
side  of  the  momentum  equations  using  the  Rossby-Haurwits. 
velocity  field  valuer  and  com pari ng  these  ayainst  analytic 
values  for  the  terms .  A  more  accurate  fourth-order  accuracy 
finite-  d.i  fforenee  scheme  will  be  devised  and  implemented  in 
the  one-layer  model  in  order  to  determine  the  role  of  the 
finite  difference  scheme  in  forecast  accuracy.  Its  values 
for  the  momentum  equation  terms  will  also  be  compared 
ayainst  their  analytic  values  usiny  Rossby -Haurwi tz  values. 
The  model  will  be  further  modified  to  accommodate  doubliny 
and  quadrupling  the  horizontal  resolution,  as  well  as  to 
include  more  complex  physical  processes  as  friction,  dia¬ 
betic  heatiny,  topoyraphy,  and  water  phase  changes.  The 
model  will  also  be  expanded  vertically  to  include  computa¬ 
tions  for  the  entire  vertical  extent  of  the  model  atmos¬ 
phere  as  shown  in  Fiy.  A-1. 
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B.  A  COMPARATIVE  STUDY  OF  WIND  SETS  FOR  AUTOMATED 
WEATHER  SYSTEMS 

I.  INTRODUCTION 

A  primary  goal  of  the  Air  Force’s  research  and 
development  program  is  to  develop  automated  weather 
observing  systems  which  will  function  unattended  in  an 
all-weather  environment.  As  a  step  toward  achieving 
this  goal,  a  study  is  underway  to  examine  four  wind  sets 
and  their  response  characteristics  under  a  variety  of 
weather  conditions. 

The  four  instruments  being  tested  are  the  Rosemount 
Orthogonal  Airspeed  System,  the  Ciimatronics  Wind  System  - 
Mark  I,  the  R.  M.  Young  Propvane,  and  the  J-Tec  Vortex 
Anemometer.  Based  on  AFGL's  extensive  prior  experience 
with,  and  confidence  in,  the  Ciimatronics  Wind  System  - 
Mark  I,  it  is  being  used  as  the  standard  of  comparison 
in  the  study. 

II.  DESCRIPTION  OF  INSTRUMENTS 

The  four  wind  instruments  were  selected  because 
each  makes  use  of  a  different  principle  for  determining 
wind  speed  or  direction.  The  J-Tec  VA-100  Vortex  Anemo¬ 
meter  consists  of  a  fiberglass  wind  vane  which  determines 
wind  direction  through  the  use  of  a  potentiometer. 

Mounted  on  the  wind  direction  vane  is  an  obstacle  which 
creates  a  series  of  vortices  in  the  downstream  flow. 

These  vortices  are  detected  by  an  ultrasonic  beam,  their 
frequency  is  determined  and  is  converted  to  wind  speed. 

The  instrument  was  designed  for  unattended  operation  in 
a  high-wind  environment. 

The  Ciimatronics  Wind  System  -  Mark  1  has  been  used 
in  a  number  of  mososcale  weather  field  studies  and  in  the 
Modular  Automated  Weather  System  (MAWS)  demonstration  at 
Scott  AFB.  It  consists  of  a  sot  of  stain  Jess  stool  onps 
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which,  in  response  to  wind,  rotate  and  cause  a  light  source 
to  be  chopped  at  a  frequency  which  is  proportional  to  the 
speed.  The  wind  direction  is  determined  by  a  vane  which 
is  connected  to  a  frictionless  direction  transducer. 

The  R.  M.  Young  Propvane  Model  8002C  consists  of  a 
wind  vane,  which  determines  wind  direction  by  use  of  a 
potentiometer,  coupled  with  a  propeller,  which  activates 
a  D.C.  generator  whose  voltage  output  is  directly  propor¬ 
tional  to  wind  speed. 

The  Rosemount  Model  853  Orthogonal  Airspeed  System 
measures  two  orthogonal  wind  speed  vectors  of  air  flow. 

The  instrument  has  no  moving  parts.  Instead,  pressure 
ports  located  on  the  sensor's  cylindrical  body  develop 
differential  pressures  which  are  proportional  to  the  wind 
speed.  Resolution  of  these  pressure  differentials  through 
a  transducer  allows  the  computation  of  wind  direction  and 
speed . 


III.  WEATHER  TEST  FACILITY 

The  test  Is  being  conducted  at  the  AFGL  Weather  Test 
Facility  at  Otis  AFE,  MA.  The  facility  was  designed  pri¬ 
marily  to  serve  as  a  test-bed  for  the  comparison  of  new 
sensors  for  Air  Force  use  and  for  the  examination  of  wea¬ 
ther  automation.  It  was  felt  that  the  wide  variety  of 
weather  that  occurs  in  this  area  would  provide  an  excellent 
environment  in  which  to  test  the  durability  and  responsive¬ 
ness  of  the  candidate  wind  sets.  Fig.  B-l  shows  the  instru¬ 
mented  towers  A,  B,  C,  P,  and  Q,  and  the  ground  site  X, 
where  the  wind  sets  are  installed. 

The  Modular  Automated  Weather  System  is  also  located 
at  the  Weather  Test  Facility  and  is  connected  to  Hanscom 
AFB  by  commercial  telephone  lines.  At  the  base  of  each 
tower  and  at  the  ground  site  (see  Fig.  B-l)  a  remote  data 
unit  (RDU)  is  located.  These  RDU’s  serve  as  links  between 
the  weather  sensors  and  the  MAWS  data  acquisition  system. 

The  RDU's  are  microprocessors  which  poll  the  sensors, 
convert  the  analog  signals  to  digital,  sum  the  values,  and 
produce  a  one-minute  mean  for  each  sensor.  Included  are 
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visibility  instruments  such  as  scatter  meters,  transmis- 


someters,  videograph,  and  nephelometer ;  and  other  instru¬ 
ments  such  as  the  wind  sets  described  in  Section  II, 
radiometers,  temperature  and  dewpoint  sensors,  rain  gauge, 
and  rotating  beam  ceilometers. 

Most  sensors  report  at  12-second  intervals;  the  wind 
sets  report  each  six  seconds.  Data  from  the  RDU's  are  fed 
to  a  supervisor  at  the  rate  of  120  characters  per  second 
for  further  processing.  Refined  means,  collected  from  all 
activated  sensors,  or  the  raw  voltages  from  up  to  20 
selected  sensors,  are  sent  to  AFGL  for  real-time  print 
display  or  archiving  on  magnetic  tape. 

A  CDC  6600  computer  program  has  been  developed  to 
process  data  from  these  magnetic  tapes  to  form  a  data 
base  for  several  research  projects,  including  cloud  base 
height  studies  and  wind  shear  investigations.  In  the  case 
of  the  wind  data,  one-minute  mean  u,  v  components  are  cal¬ 
culated,  from  which  the  resultant  one-minute  mean  wind 
direction  and  speed  are  obtained.  The  one  exception  to 
this  procedure  are  the  data  from  the  Rosemount  Orthogonal 
Airspeed  System,  which  are  initially  obtained  in  u,  v 
component  form. 

IV.  DATA  ANALYSIS 

Software  has  been  developed  to  analyze  the  data 
collected  by  the  four  wind  sets.  Two  correlation  pro¬ 
grams  were  written.  The  first,  which  is  executed  as  part 
of  the  routine  analysis  of  the  Otis  data,  compares  the 
wind  direction  and  wind  speed  from  the  wind  sets  by  means 
of  scatter  diagrams.  In  the  comparison  of  the  Rosemount 
with  the  Climatronics ,  three  wind  speed  regimes  (wind 
speed  greater  than  4.0  m/sec,  greater  than  2 .  'j  m/sec,  and 
all  wind  speeds)  are  evaluated  due  to  known  characteristics 
of  the  Rosemount.  The  scatter  program  also  compares  the 
J-Tcc  and  the  R .  M.  Young  with  the  Climatronics,,  using  the 
all-wind -speed  category.  In  adc.it ion  to  1  he  computer- 
printed  scatter  plots,  calculations  arc  made  which  riv< 
the  equation  of  best  fit,  the  HMD  error,  the  corre 1  at i on , 


an  average  of  each  variable,  the  standard  deviations, 
the  co-variance,  and  the  standard  error  of  estimate. 

The  second  program  produces  a  pen  and  ink  plot  of 
the  scatter  diagrams,  as  illustrated  in  Figs.  B-2 
through  B-5.  Information  displayed  on  the  plots  consists 
of  the  date-time  group  (e.g.,  6  September  1979).  "the  num¬ 
ber  of  observations  (N),  and  the  equation  of  the  line  of 
best  fit.  These  figures  are  provided  to  illustrate  how 
the  scatter-  diagram  analysis  is  being  used  in  the  evalua¬ 
tion  of  the  four  wind  sensors.  Fig.  B-2  shows  that  the 
J-Tec  wind  speed  is  consistently  higher  than  the  Clima- 
tronics;  therefore,  recalibration  of  the  voltage- to-wind 
speed  algorithm  would  be  required.  However,  the  relation¬ 
ship  is  linear  and  there  is  very  little  scatter.  The  wind 
direction  comparison  {Fig.  B-3)  covers  a  small  range  of 
direction  but  the  fit  is  excellent.  The  wind  speed  com¬ 
parison  between  the  Climatronics  and  the  R.  M.  Young 
(Gill)  (Fig.  B-4)  shows  that  the  R.  M.  Young  consistently 
reports  lower  speeds  than  the  Climatronics,  and  again 
algorithm  adjustment  is  in  order,  but  the  relationship 
is  linear  and  the  scatter  is  quite  small.  The  wind  direc¬ 
tion  comparison  between  the  two  (Fig.  B-5)  is  very  good. 

In  the  later  phases  of  the  study,  after  a  sufficiently 
large  data  set  has  been  accumulated,  the  calibration 
algorithms  will  be  adjusted  prior  to  the  comparison  based 
on  weather  conditions. 

In  order  to  compare  the  wind  sets  in  different  weather 
conditions,  use  will  also  be  made  of  a  decision- tree  com¬ 
puter  program  developed  for  another  project.  The  decision- 
tree  program,  using  data  from  the  automated  instrument 
array  at  Otis,  objectively  determines  the  kinds  of  weather 
and  obstructions  to  vision  that  are  occurring.  Fig.  B-6 
lists  the  elements  normally  observed  by  human  observers 
and  indicates  those  which  the  decision-tree  program  has, 
the  ability  to  specify.  This  program  is  being  used  to 
classify  the  periods  during  which  the  wind  sets  will  be 
compared . 
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C.  Me I DAS  HARDWARE  AND  SUPPORT  TO  RESEARCH  PROJECTS 

I.  INTRODUCTION 

SASC‘s  objective  is  to  improve  the  output  of  the 
AFGL/Meteorology  Division  McIDAS  Facility  and  to  operate 
it  in  support  of  meteorologi cal  research.  As  an  inter¬ 
active  system  McIDAS  combines  its  resident  programs  and 
processing  speed  with  the  skill  and  ingenuity  of  the  user, 
to  facilitate  handling  large  masses  of  data.  It  acquire;', 
processes,  and  analyzes  satellite  and  conventional  weather 
data.  Visible  and  infrared  satellite  data  are  received 
from  two  geosynchronous  satellites  operated  by  NOAA/NESS. 
Conventional  data  (surface  and  upper  air  reports)  arrive 
via  high  speed  FA A  data  line. 

II.  Me  I  DAD  HARDWARE  ENHANCEMENT 

MeTDAS  consists  of  a  real-time  data  acquis i tion  sub¬ 
system,  Harris  r Qt’A/e  computer  with  PSK  memory  (fA-bit 
word  length),  80  megabyte  digital  dirk  drive,  four  addi¬ 
tional  d i sk  drives  each  with  10  megabyte  capacity,  line 
printer,  card  render,  two  tape  transports  for  track  800 
rid  magnetic  tapes,  1 ! C  video  disk  recorder  with  a  total 
ra  par  •  !y  of  tPO  t’rames,  and  two  t'KT,  keyboard  eonso  j  es  for 
the  sys  tern  of  ■■■  ra  t  err. . 

A  mnjrr  system  im>  rovriii'-td  wars  1  lie  i  r.s  ta  i  la  ‘  :  o!  r  !' 


!  he  s’  hyt.e  digital  disk,  whiej  !:ct  on  ly  l  rev  . 

•id  !  >  t  ’em  !  .  ■  f  e  e;  i  .  .  \a  l  :  I  e  1  t  ",  hit  can  a  e ,  a  s  S  da  ‘a  a  * 


as  transmitted  every  half-hour  from  the  satellite.  Sony 
archive  recorder  components  include  a  modified  VO-285OA 
video  cassette  recorder,  a  motor  drive  amplifier  to  run 
the  recorder  at  500  rpm,  and  the  encoding  electronics  to 
format  GOES/SMS  data  for  recording  on  tape.  Player  compo¬ 
nents  are  a  modified  VF-1200  video  cassette  player,  a 
motor  drive  amplifier  to  run  the  player  scanner  at  500 
rpm,  an  Aydin  Monitor  Model  ^50  bit  synchronizer,  and 
decoding  electronics  to  reconstruct  GOES/SMS  data  as  well 
as  to  provide  search  information  and  remote  control  capa¬ 
bility.  The  University  of  Wisconsin  Space  Science  and 
Engineering  Center  assembled  the  system,  to  which  a 
switching  matrix  and  byte  mangier  were  added  to  integrate 
fully  the  archive  units-  with  Me  IDAS .  Two  archive  systems 
were  installed;  one  on-line,  the  other  as  a  spare  in  anti¬ 
cipation  of  simultaneous  recording  from  two  satellites. 

MeTDAG  data  lines  were  removed  and  replaced  with  a 
full  duplex  dataphor.e  system.  The  Dataset  212A  is  capable 
of  t  ransm  i  1 1.  i  ng  and  receiving  data  at  a  rate  of  1,200  BPS. 
System  compatibility  with  the  University  of  Wisconsin  in¬ 
stallation  now  makes  ox change  of  software  technically 
fear  it  le. 

Modi  float  ions  were  made  in  the  cursor  control 
electronics  in  an  effort  to  stabilize  the  cursor  position 
generated  by  tic  „oyst  i  cks .  fa-design  employed  a  +18  volt 
and  -8  volt  isolated  power  suj t iy  instead  of  the  +18  volts 
swrplied  from  the  system.  The  isolated  power  supply , 
liens  with  resistors  on  b<  !  h  •  t.ds  of  the  joysticks-,  do- 
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III.  McIDAS  SUPPORT  TO  RESEARCH  PROJECTS 


McIDAS  supported  a  series  of  ABRES  (Advanced  Ballistic 
Re-entry  System)  missions  conducted  at  the  Kwajalein  Missile 
Range  (KMR).  Missions  were  weather-dependent,  categorized 
as  heavy  weather,  moderate  weather,  or  clear  air,  accord¬ 
ing  to  the  conditions  desired  for  re-entry. 

Because  of  the  isolation  of  Kwajalein  and  the  absence 
of  surrounding  surface  stations,  local  weather  observers 
depended  heavily  on  radar  and  ship/aircraft  reports  when 
making  forecasts.  Prior  to  mission  time,  McIDAS  provided 
KMR  forecasters  with  a  single  image  showing  cloud  top 
temperatures,  a  cross-correlation  of  two  or  more  images 
that  yielded  wind  vectors,  and  loops  of  several  images 
that  depicted  cloud  growth  and  decay. 

Another  valuable  input  to  ABRES  was  the  Environmental 
Severity  Index  (ESI),  which  was  a  measure  of  particulate 
matter  at  various  levels  in  the  atmosphere.  The  measure¬ 
ment  indicated  the  potentially  erosive  effects  of  such 
matter  on  a  re-entry  vehicle.  On  several  occasions 
McIDAS  was  instrumental  in  determining  whether  a  mission 
count  down  was  to  continue  to  launch  or  not. 

The  Large  Scale  Cloud  System  Program,  a  continuing 
AFGL  effort  to  obtain  better  definition  of  the  micro¬ 
physics  of  large  winter  storms,  sampled  extensive  homo¬ 
geneous  cloud  masses  as  they  progressed  from  west  to  east 
over  North  America.  Aircraft  data  relating  to  particle 
size,  type,  number,  and  liquid  water  content  were  collected 
at  various  atmospheric  levels.  McIDAS  satellite  imagery 
loops  showing  cloud  motion  and  cloud  top  temperatures 
enabled  the  ground  support  group  to  identify  the  most 
desirable  sampling  locations  to  the  inflight  crew. 

Because  of  the  high  cost  of  flight  operations  it  was  impor¬ 
tant  for  the  aircraft  to  reach  the  best  location  as. 
quickly  as  possible.  McIDAS  provided  conventional  surface 
observation  plots  and  analyses  of  wind  vectors,  pressure- 
fields,  streamlines,  dew  point  temperatures,  precipitation 
amounts,  current  weather,  visibility,  and  ciond  (-over  con¬ 
ditions  (  see  Kirs.  C-P ,  f-'i,  C-  h  ,  arid  h’-H  ) .  A  Iso  ,  up|  or 
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Figure  C-2a.  All  Ohio  surface  station  reports;  (lower) 
list  of  individually  selected  stations. 
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Figure  C-Ob.  All  surface  stations  n-iori  inr  ra  i : 


Figure  C-4a.  Surface  pressure  analysis  for  the  U.S.  The 
20  isobar  =  1,020  mb;  -4  isobar  =  996  mb. 


Figure  C--  b.  Surface  analysis  of  wind  flow  (streamlines) 
centered  on  Minnesota. 


Figure 


Figure 


-5a.  500  mb  station  model  plot  centered  over  Ohio. 


-5b.  300  mb  wind  vectors  for  the  United  States. 
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air  height  analyses,  vorticity  analyses,  and  radiosonde 
(Stuve)  and  cross-section  plots  (see  Pigs.  C-5,  C-6,  C-7) 
helped  determine  suitable  sampling  altitudes.  Data  were 
then  archived  on  tape  and  hard-copy  printouts  for  later 
documentation  of  the  synoptic  situation. 

Another  on-going  study  requiring  McIDAS  support'  was 
the  Mesoscale  Forecasting  Program.  It  was  concerned  with 
developing  short  range  forecast  techniques  using  sequences 
of  GOES/SMS  imagery.  A  microprocessor-controlled  tape 
drive  known  as  ODIS  (Off-line  Data  Ingestion  System;  see 
Fig.  C— 1 )  collected  one-mile  visible,  four-mile  infrared 
images  centered  over  eastern  Pennsylvania,  at  hourly 
intervals  from  14002  to  2100Z.  Data  were  collected  at 
half-hourly  intervals  if  the  study  area  contained  parti¬ 
cularly  interesting  or  severe  weather.  Collected  data 
were  then  saved  on  magnetic  tape,  three  days'  data  per 
tape.  The  data  were  analyzed  on  the  AFGL  central  com¬ 
puter  to  yield  motion  vectors,  brightness,  and  IR  tempera¬ 
tures  which  were  converted  into  weather  parameters  for 
the  forecast  procedures. 

The  AFGL  Icing  Program  was  designed  to  observe  and 
measure  accumulation  of  ice  on  an  aircraft.  Five  flights 
were  conducted  with  an  instrumented  C-I30  aircraft, 
which  sampled  over  any  one  of  12  radiosonde  stations 
within  a  300-mile  radius  of  Wright-Patterson  AFB,  Ohio. 
Sequences  of  visible  and  IR  imagery  were  provided  by 
McIDAS  to  determine  if  suitable  icing  conditions  existed 
within  the  study  area.  Other  data  included  hourly  FAA 
surface  observations  of  current  weather  (see  Fig.  C-8), 
temperature,  dew  point  temperature,  pressure,  wind 
speed  and  direction,  cloud  cover  and  ceiling,  6-hour 
precipitation  total,  streamlines  (see  Figs.  C-2,  C-3, 

C-4),  upper  air  height  analyses,  radionsonde  and  cross- 
section  plots  (see  Figs.  C-5,  C-6,  C-?).  Again,  pre- 
flight  and  in-flight  data  were  archived  on  tape  and 
printouts  for  later  inspection  at  report  time. 


IV.  DATA  NAVIGATION  IMPROVEMENT 


When  working  with  satellite  data  it  is  important  to 
know  precisely  the  latitude  and  longitude  of  any  given 
point  within  the  field  of  view.  The  situation  is  exacer¬ 
bated  if  the  imagery  contains  an  expanse  of  ocean  or  cloud 
obscured  land.  To  insure  maximum  location  accuracy  of 
satellite  data,  a  navigation  program  was  developed  at  the 
University  of  Wisconsin  Space  Science  and  Engineering 
Center.  That  program  relates  a  simplified  orbital  model 
in  McIDAS  to  the  scan  line  and  element  of  the  satellite. 
Orbital  data  used  in  the  navigation  process  are  computed 
by  NOAA  from  a  complex  model  taking  into  account  several 
celestial  bodies.  Since  the  McIDAS  orbital  model  is  only 
a  two-body  system  it  frequently  requires  adjustment  when 
a  manuever  (re-orientation)  occurs  or  when  an  orbit  decays 

Real-time  data  used  in  the  navigation  procedure  (see 
Fig.  C-9)  are  0.5  mi  resolution  images,  enlarged  on  the 
TV  screen  by  a  factor  of  four.  The  landmark  is  generally 
a  prominent  point  along  a  coastline  with  a  distinct  land/ 
sea  interface.  After  five  to  seven  images  have  been  in¬ 
gested,  the  landmark  common  to  all  is  defined;  that  is, 
the  latitude  and  longitude  of  the  particular  point  are 
called  up  from  a  locate  file  and  entered  into  the  system. 
An  acetate  overlay  map  is  placed  over  the  TV  screen,  with 
a  cross-hair  ^entered  on  the  exact  point.  Once  image  and 
overlay  are  aligned,  a  digital  cursor  on  the  TV  screen  is 
matched  to  the  cross-hair  and  the  satellite  scan  line  and 
element  are  measured.  Values,  from  each  image  are  com¬ 
pared  with  the  line  and  element  numbers  derived  from  the 
Me T DAD  model.  Differences  between  the  two  sets  are  known 
as  the  residuals.  If  the  residuals  are  small  (less  than 
±  ? .  c,0 ,  0  being  perfect)  it  is  considered  a  good  orbit. 
Large  residuals  indicate  that  adjustment  of  one.  or  more  of 
the  NOAA  orbital  parameters  is  necessary.  Extremely  large 
residuals  indicate  that  a  new  orbit  may  he  necessary . 

Tri  1 979  LA  DC  initiated  a  project  to  develop  a  logical 
methodology  for  adjusting  orbital  parameters. 

As,  a  brief  example,  Fig.  C-10  is.  a  navigation  print- 
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Fifiirp  C-l  0 .  Navigation  printout  from  Novombor  28,  1979. 
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out  from  November  28,  1979*  From  the  residuals  it  is 
apparent  that  the  orbit  for  that  day  was  discrepant. 

Minor  adjustment  was  sufficient  to  bring  the  residuals 
within  acceptable  range.  One  orbital  parameter  whose  in¬ 
crease  or  decrease  greatly  affects  the  residuals  is  the 
semi-major  axis,  which  is  defined  as  one-half  the  greater 
axis  of  the  ellipse  of  motion.  The  semi-major  axis  in 
Fig.  C-10  is  42,167.62  km.  Reducing  it  by  3  km  to 
42,164.62  (see  Fig.  C-ll)  reduces  the  residuals  to 
acceptable  values.  Thus,  given  any  point  within  the 
landmark  images,  its  location  would  be  less  than  one 
nautical  mile  in  error.  Although  adjusting  the  orbital 
parameters  is  not  done  daily,  the  imperfect  satellite 
orbit  requires  that  navigation  be  done  daily.  A  good 
orbit  should  last  approximately  two  weeks  before  the 
residuals  become  excessive. 

Other  orbital  elements  -  eccentricity,  inclination, 
mean  anomaly,  argument  of  the  perigee,  and  right  ascen¬ 
sion  -  remain  to  be  investigated. 
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D.  McIDAS  SYSTEMS  AND  APPLICATIONS  SOFTWARE 


I.  INTRODUCTION 

McIDAS  software  improvement  involved  two  types  of 
activities  -  modifying  existing  programs  and  adding  new 
ones.  The  first  included  correcting  logic  faults  and  key¬ 
punching  errors,  increasing  utility  of  individual  pro¬ 
grams  by  adding  new  functions,  and  adjusting  programs  to 
accommodate  extensive  changes  in  the  hardware  system. 
Certain  new  programs  were  supplied  by  the  University  of 
Wisconsin  and  adapted  for  the  AFGL  McIDAS,  while  others 
were  written  by  SASC  to  meet  specific  needs  of  AFGL  and 
SASC  scientists. 

II.  MODIFICATIONS 

By  early  1979  all  system  software  and  most  applica¬ 
tions  software  originally  provided  by  the  University  of 
Wisconsin  had  been  implemented.  However,  one  category  of 
programs  called  NOWCASTING  still  contained  logic  errors. 
Since  these  programs,  designed  to  analyze  and  display- 
conventional  weather  data,  are  used  routinely  for  monitor¬ 
ing  current  weather,  such  faults  quickly  became  apparent 
and  were  corrected.  Additional  functions  were  included 
to  expedite  location  and  identification  of  individual 
surface  and  radiosonde  stations,  and  a  plot  file'  capa¬ 
bility  was  added  to  those  programs  written  for  WRRRM 
(Write  Random  Read  Raster  Memory)  display  only.  This 
modification  makes  it  possible  to  write  into  the  display 
file  and  transfer  the  finished  product  to  a  TV  frame. 

The  mod i fi ca t ions  required  to  make  the  programs 
written  for  the  University  of  Wisconsin  McIDAS  compatible 
with  the  AFGL  system  were  due  primarily  to  cor  'igurational 
differen-es  between  the  two  facilities.  One  major  differ¬ 
ence,  for  example,  involved  the  form  in  which  support 
software  was  stored  on  disk.  Most  AFG1  support  routines 
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existed  as  individual  relocatable  modules,  whereas  the 
University  of  Wisconsin  had  gathered  them  into  special¬ 
ized  libraries.  Not  only  does  the  latter  procedure  reduce 
the  amount  of  disk  storage  required  but  also  simplifies 
the  job  control  stream  of  most  applications  programs. 

AFGL  software  was  therefore  reorganized  to  incorporate 
this  improved  structure. 

The  addition  of  the  Sony  video  cassette  archive  ne¬ 
cessitated  adaptation  of  the  real-time  satellite  ingest 
program  and  support  routines  to  permit  ingestion  of 
imagery  from  the  archive  playback  unit.  Two  test  programs 
provided  by  the  University  of  Wisconsin  for  sampling  the 
quality  of  recorded  imagery  were  also  implemented. 

Increase  in  mass  storage  afforded  by  the  addition  of 
an  80  megabyte  disk  made  feasible  an  increase  from  16  to 
24  in  digital  areas  available  to  each  user  terminal  for 
storage  of  satellite  imagery.  This  increase  required 
creation  of  new  data  files  and  adaptation  of  every  appli¬ 
cations  program,  system  routine,  and  library  module  which 
must  have  access  to  them. 

III.  NEW  PROGRAMS 

Additional  storage  space  also  became  available  for 
new  programs.  The  Area  Statistics  software  package,  from 
the  University  of  Wisconsin,  is  a  collection  of  programs 
and  data  files  designed  to  select  portions  of  satellite 
images  and  perform  statistical  analyses  on  the  digital 
data  they  represent.  Although  these  programs  are  still 
in  the  developmental  stage,  functions  which  are  available 
include  calculation  of  brightness  distributions  in  histo¬ 
gram  form,  surface  areas,  and  area  sums  lying  between, 
above,  below,  or  at  specific  brightness  levels. 

The  largest  addition  to  McIDAS  software  was  a  collec¬ 
tion  of  relocatables ,  data  files,  and  a  main  driver  pro¬ 
gram,  designated  the  Mesoscale  Forecasting  Facility  (MFF). 
Motivation  for  this  major  programming  effort  was  the  Air 
Force  desire  to  make  the  many  valuable  McIDAS  forecasting 
aids  available  to  all  AFGL  Mesoscale  Forecasting  Branch 
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meteorologists  and  to  encourage  regular  exercise  of  short- 
range  prediction  skills. 

IV.  MESOSCALE  FORECASTING  FACILITY  (MFF) 

With  only  minimal  knowledge  of  McIDAS  operating  pro¬ 
cedures,  a  forecaster  using  MFF  can  access  most  of  the  NOW- 
CASTING  programs  and  many  satellite  imagery  request  and  dis¬ 
play  functions,  formulate  a  forecast,  enter  his  predictions 
into  a  personalized  data  file,  verify  the  forecast  at  a 
later  time,  and  evaluate  his  performance  by  reviewing  his 
past  forecasts  and  accumulated  skill  scores. 

MFF  is  a  menu-driven  program;  that  is,  instruction 
pages  present  a  series  of  multiple-choice  questions. 

Answers  provided  by  the  forecaster  are  used  to  build  normal 
McIDAS  commands  performing  a  variety  of  functions  which 
assist  in  formulating,  inputting,  and  verifying  short-range 
(1-6  hr)  forecasts. 

The  MFF  program  consists  of  four  major  components  - 
chief  forecaster,  operational  forecaster,  interrogation, 
and  verification  modules.  The  chief  forecaster  module  en¬ 
ables  the  responsible  scientist  to  select  the  station  for 
which  the  day's  forecasts  are  to  be  made  and  the  parameters 
(predictands)  whose  values  are  to  be  forecast.  He  also  leads 
a  brief  discussion  of  current  weather  conditions.  Opera¬ 
tional  forecasters  participate  in  the  weather  discussion  and, 
using  the  operational  forecaster  module,  formulate  fore¬ 
casts  in  compliance  with  the  specifications  set  by  the 
chief.  The  interrogation  module  permits  the  operator  to 
review  his  performance  to  date  by  listing  past  forecasts 
and  accumulated  skill  scores.  The  verification  module  com¬ 
pares  the  predictions  of  each  forecaster  with  the  actual 
station  reports  and  computes  individual  skill  scores. 

A.  CHIEF  FORECASTER  MODULE 

The  primary  component  of  MFF  is  the  module  for 
the  chief  forecaster.  After  careful  examination  of  the  most 
recent  weather  data,  including  conventional  surface  and  upper 
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air  reports  and  real-time  satellite  imagery,  the  chief  fore¬ 
caster  chooses  a  region  exhibiting  interesting  meteorologi¬ 
cal  features  and  provides  a  selection  of  challenging  fore¬ 
cast  parameters. 

Three  media  are  available  for  conventional  data  repre¬ 
sentation.  The  most  flexible  method  utilizes  the  video 
graphics  display  capabilities  of  McIBAS.  Plots  and  con¬ 
toured  analyses  of  surface  and  upper  air  data  are  produced 
in  three  colors  on  a  TV  monitor.  Values  of  any  reported 
parameters  and  some  additional  derived  parameters  (Tabs. 

D-l  and  D-2)  may  be  plotted  over  U.S.  or  regional  base  maps. 
Fig.  C-3a  is  an  example  of  a  surface  temperature  plot  over 
New  England.  In  addition  to  digital  data,  special  symbols 
are  available  for  representation  of  sky  cover  and  weather 
(Fig.  C-3b).  Wind  speed  and  direction  can  be  displayed  in 
three  forms s  digitally  as  in  the  station  model  plot  depicted 
in  Fig.  C-5a,  vectorially  as  in  Fig.  C-5b,  and  as  standard 
meteorological  flags.  Those  parameters  which  exist  as  con¬ 
tinuous  or  nearly  continuous  fields  (Tabs.  D-3  and  D-4) 
can  be  interpolated  to  grid  point  values  and  objectively 
contoured.  Figs.  C-4  and  C-6  show  examples  of  surface  and 
upper  air  analyses  performed  In  this  manner  and  displayed 
by  McIDAS,  Two  additional  types  of  analyses  are  available 
for  radiosonde  reports.  Individual  station  reports  in  the 
form  of  Stuve  thermodynamic  diagrams  (Fig.  C-7a)  are  useful 
in  estimating  stability.  Groups  of  four  to  six  reports 
organized  in  spatial  or  temporal  series,  analyzed  as  verti¬ 
cal  c^oss-sections  or  time  sections,  provide  a  detailed 
view  of  the  distribution  of  thermal,  moisture,  and  velocity 
fields  in  the  atmosphere  (Fig.  C-7b). 

The  time  required  for  McIDAS  to  generate  any  of  these 
types  of  displays  is  normally  less  than  two  minutes,  de¬ 
pending  on  the  extent  to  which  the  system  is  involved  in 
other  activities  concurrently.  As  leader  of  the  weather 
discussion,  the  chief  may  find  that  this  time  delay 
severely  limits  the  number  of  parameters  he  is  able  to 
present.  His  module  therefore  provides  the  option  to  save 
the  specifications  for  up  to  five  commands  generated  while 
preparing  his  presentation.  A  minor  amount  of  time  is 


Table  D -1 .  Parameters  Available  for  Surface  Data  Plots 


Temperature 

Dew  point  temperature 

Pressure 

Wind  direction  and  speed 
Low  cloud  cover 
Middle  cloud  cover 
High  cloud  cover 
Current  weather 
♦Current  weather  symbols 


♦Wind  flags 
Visibility 
Wind  gusts 
Precipitation 
Low  cloud  height 
Middle  cloud  height 
High  cloud  height 
♦Ceiling  height 
♦♦Potential  temperature 
♦♦Equivalent  potential 
temperature 


♦Not  available  on  line  printer  plots 
♦♦Available  on  line  printer  only 


Table  D-2 .  Parameters  Available  for  Upper  Air  Data  Plots 


Station  identifiers 
Temperature 
Dew  point  temperature 
Wind  direction  and  speed 
Heights 
♦Wind  flags 
♦Wind  vectors 

Plot  all  parameters  in  station  model  format 


♦Not  available  on  line  printer  plots 


Table  D-3*  Parameters  Available  for  Surface  Data  Contouring 


Temperature 

Dew  point  temperature 

Pressure 

Potential  temperature 
Equivalent  potential 
temperature 
Mixing  ratio 
Wind  speed 
Total  cloud  cover 
Divergence 
Vorticity 

Temperature  advection 


Dew  point  advection 
Mixing  ratio  advection 
Streamlines 
Ceiling  height 
Stretching  deformation 
Shear  deformation 
Dew  point  divergence 
Mixing  ratio  divergence 
Visibility 
Low  cloud  cover 
Middle  cloud  cover 
High  cloud  cover 
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Table  D-4.  Parameters  Available  for  Upper  Air  Data  Contouring 

Temperature  Vorticity 

Dew  point  temperature  Temperature  advection 

Heights  Dew  point  advection 

Potential  temperature  Mixing  ratio  advection 

Equivalent  potential  Streamlines 

temperature  Stretching  deformation 

Mixing  ratio  Shear  deformation 

Wind  speed  Dew  point  divergence 

Total  cloud  amount  Mixing  ratio  divergence 

Divergence 

still  required  to  retrieve  and  analyze  the  data  and  draw 
the  display  on  the  TV  screen.  If  more  than  five  charts  are 
desired,  or  if  the  time  allotted  for  producing  each  one 
must  be  further  reduced,  a  second  display  medium  is  avail¬ 
able.  Any  of  the  above  types  of  displays  may  be  written 
into  a  plot  display  file  and  transferred  to  the  analog 
disk  where  they  are  preserved  for  future  use.  Images  pro¬ 
duced  in  this  manner  forfeit  the  advantage  of  color. 

The  third  display  medium  provides  a  means  of  producing 
hard  copies  of  most  of  the  charts  described.  Line  printer 
output  is  requested  in  lieu  of  TV  display.  Surface  and 
upper  air  data  plots  are  printed  over  a  U.S.  base  map. 

Fig.  D-l  shows  a  line  printer  plot  of  current  weather. 

(Special  weather  symbols  seen  in  Fig.  C-3b  are  not  avail¬ 
able  through  this  device.)  The  temperature  plot  in  Fig. 

D-2  illustrates  the  approximate  coverage  of  surface  data. 

Upper  air  coverage  is  shown  in  Fig.  D-3  by  a  plot  of  the 
five-digit  identifiers  of  all  reporting  radiosonde  stations. 
Contouring  of  data  is  simulated  on  line  printer  output  by 
alternate  shading  between  contours.  Examples  of  surface 
analyses  of  this  type  are  shown  covering  the  U.S.  in  Fig. 

D-4  and  on  a  regional  scale  in  Fig.  D-5.  Upper  air  charts 
analyzed  in  this  fashion  are  also  available.  Examples  of 
500  mb  heights  and  vorticity  are  given  in  Figs.  D-6  and 
D-7  respectively.  While  soundings,  cross-sections,  and  time 
sections  cannot  be  provided  in  analyzed  form  by  the  line 
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F ini  re  D-l .  Line  printer  plot  of  current  weather. 
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rir.ter  plot  of  radiosonde  station  identifiers. 
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igure  D-5.  line  printer  analysis  of  temperature  over  New  England 
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printer,  data  can  be  listed  to  facilitate  hand  plotting 
(Fig.  D-8)  or  plotted  for  hand  analysis  (Fig.  D-9). 

If  the  operator  indicates  a  desire  to  examine  conven¬ 
tional  data,  MFF  will  lead  him  through  a  series  of  questions 
and  answers  until  enough  information  has  been  elicited  for  a 
display  to  be  generated.  The  operator  first  chooses  the 
output  medium  most  suited  to  his  needs.  Then  a  choice  is 
made  of  type  of  data  (surface  or  upper  air)  in  which  he  is 
interested  and  the  form  in  which  it  is  t o  be  displayed.  A 
specific  parameter  must  be  selected  for  all  but  those  dis¬ 
plays  representing  the  entire  report  for  a  given  station. 

The  time  period  during  which  the  data  were  taken,  hourly 
for  surface  reports  and  twice  daily  (0000  and  1200  GMT) 
for  upper  air  soundings,  is  then  specified.  For  radiosonde 
data  a  quasi-horizontal  pressure  surface  for  which  the 
analysis  is  to  be  performed  is  also  specified;  it  must  be 
one  of  the  mandatory  reporting  levels.  The  last  major 
option  is  the  map  projection  as  background  for  the  analy¬ 
sis.  Choices  range  from  a  full  U.S.  map  to  regional  cover¬ 
ages  centered  on  any  of  the  50  states.  Several  additional 
options  are  provided  for  determining  color,  contour  inter¬ 
val,  and  display  format. 

In  addition  to  conventional  data  the  chief  forecaster 
may  request  ingestion  of  imagery  from  one  of  the  two  geo¬ 
stationary  satellites,  SMS-2  located  off  the  East  coast  at 
75°  West  longitude  and  G0ES-3  off  the  West  coast  at  135°  W. 
Each  satellite  transmits  full  disk  imagery  at  half-hour  in¬ 
tervals.  SMS-2  image  times  occur  on  the  hour  and  half- 
hour,  while  G0ES-3  begins  transmission  at  15  and  45  minutes 
past  the  hour.  In  submitting  imagery  requests  the  chief 
forecaster  specifies  the  region  of  interest,  preferred 
resolution,  and  image  times.  Visible  and  IR  data  are 
transmitted  simultaneously,  so  for  every  image  transmission 
ingested  a  pair  of  pictures  (one  from  each  channel)  is  pro¬ 
duced,  coinciding  in  spatial  and  temporal  coverage.  A 
series  of  up  to  three  picture  pairs  may  be  requested  for 
this  application.  The  visible  sequence  is  loaded  on  con¬ 
secutive  TV  frames  and  the  IR  sequence  is  loaded  on  the  oppo¬ 
site  channel  of  the  analog  dirk  on  consecutive  frames. 
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Images  may  then  be  examined  individually  or  as  short  movie 
loops.  It  is  also  possible  to  study  satellite  and  conven¬ 
tional  surface  data  simultaneously  by  overlaying  analyzed 
surface  reports  on  the  imagery  (Fig.  C-8). 

After  completing  his  examination  of  current  weather 
information  the  chief  determines  the  forecast  assignment  for 
the  day.  He  selects  a  surface  reporting  station  as  the  fore¬ 
cast  site  and,  from  the  list  given  in  Tab.  D~5 ,  three  pre- 
dictands.  He  also  enters  the  "persistence  forecast."  He 
initiates  the  process  with  a  request,  and  the  actual  input 
of  data  is  performed  internally  within  the  program.  This 
function  requires  retrieval  of  the  most  recent  hourly  sur¬ 
face  report  submitted  by  the  forecast  station  and  input  of 
the  reported  values  of  the  three  predictands.  Persistence 
assumes  that  these  values  will  continue  to  represent  actual 
conditions  throughout  the  forecast  period.  Persistence 
will  later  be  used  as  one  standard  against  which  the  perfor¬ 
mance  of  the  forecasters  will  be  judged. 

Finally,  the  chief  enters  his  own  forecast.  Three 
predicted  values  are  entered  for  each  predictand.  The 
first  is  a  one-hour  forecast,  the  second  verifies  in  three 
hours,  and  the  third  in  six.  Multiple  forecast  lengths 
were  selected  to  cover  a  broad  range  of  mesoscale  weather 
phenomena. 

B.  OPERATIONAL  FORECASTER  MODULE 

The  primary  function  of  the  operational  forecaster 
module  is  to  allow  the  participating  forecasters  to  enter 
their  forecasts.  Capability  to  display  conventional  weather 
data  and  satellite  imagery  is  also  provided. 

All  of  the  conventional  data  plotting,  listing,  and 
analysis  options  available  to  the  chief  forecaster  are  pro¬ 
vided  to  the  operational  forecaster.  This  enables  him  to 
check  independently  those  parameters  he  considers  necessary 
to  formulate  a  forecast  but  which  were  not  deemed  relevant 
by  the  chief  in  the  map  discussion.  The  operational  fore¬ 
caster  does  not  have  access  to  the  plot  file  medium,  saved 
commands,  or  saved  displays,  since  their  purpose  i s  to  pro- 
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vide  easy  display  of  the  plot  or  analysis  during  the  chief's 
discussion. 

The  operational  forecaster  can  display  the  satellite 
imagery  which  the  chief  requested  in  any  manner  provided  by 
the  facility;  i.e.,  individual  frames  or  visible  or  IR 
movie  loops.  He  may  also  overlay  a  data  plot  or  surface 
analysis  on  the  imagery.  This  capability  is  limited  to 
existing  imagery  since  no  imagery  request  function  is  pro¬ 
vided  by  this  module. 

The  forecast  entry  routine,  also  provided  in  the  chief 
forecaster  module,  comprises  the  major  portion  of  the  opera¬ 
tional  forecaster  module.  One  option  of  this  routine 
allows  the  forecaster  to  list  the  forecast  station  and  the 
three  predictands  being  used  that  day.  A  second  option 
initiates  the  man-machine  dialog  through  which  the  forecast 
is  entered.  Fig.  D-lOa  gives  an  example  of  the  first  option. 
The  station  is  identified  by  its  three  letter  Service-A 
identifier  and  the  predictands  are  listed  in  the  order  in 
which  they  appear  in  the  forecast  entry  dialog. 

The  second  option,  the  forecast  entry  dialog,  consists 
of  three  parts:  listing  of  the  last  seven  surface  reports 
for  the  forecast  station,  sample  forecast  entry,  ar.d  actual 
forecast  questions. 

The  first  part  of  the  dialog  displays  the  last  seven 
hourly  surface  reports  of  the  forecast  station  in  chrono¬ 
logical  order  to  help  the  forecaster  identify  significant 
trends  which  will  influence  his  forecast.  The  forecaster  is 
given  the  option  of  having  a  hard  copy  of  these  data  should 
he  desire  to  refer  to  them  later  in  the  dialog.  While  this 
display  is  being  prepared  the  forecast  entry  routine  is 
calculating  the  time  at  which  the  forecast  is  being  entered 
and  the  value  of  each  predictand  at  that  time.  These  data 
along  with  the  forecaster's  ID  arc  then  placed  into  a  buffer 
that  will  eventually  contain  the  entire  forecast. 

The  second  part  of  the  dialog  is  a  sample  forecast 
entry  which  familiarizes  the  forecaster  with  the  data  entry 
format  (see  Fig.  D-lOb).  Forecast,  data  entered  are  a  s<  t  of 
percent  probabilities  of  occurrence  t.h  <  >\  are  assigned  to 
categories  which  are  ranges,  of  values  describing  the  pre- 
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"Status  of  Forecast  Specifications"  page  from 
the  Mesoscale  Forecasting  Facility  (MFF). 
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Figure  D-lOb.  Sample  forecast  entry  page  from  MFF. 


dictand.  The  lone  exception  to  this  data  format  occurs 
when  predictand  #15,  weather  type,  is  chosen.  The  forecast 
entry  then  consists  of  a  number  associated  with  a  weather 
type  and  percent  probability  of  its  occurrence.  (See  Tab. 

D-5  for  a  listing  of  categories  associated  with  each  predic¬ 
tand,  and  Fig.  D-llb  for  a  list  of  weather  types.) 

The  sample  message  begins  by  describing  the  restrictions 
placed  on  this  set  of  probabilities  to  allow  calculation  of 
meaningful  skill  scores.  These  restrictions  are:  probabili¬ 
ties  must  be  integers  in  the  range  0  ^  x  ^  100,  a  single 
category  must  have  the  highest  probability  assigned  to  it, 
and  the  sum  of  the  probabilities  must  equal  100. 

The  remainder  of  the  message  displays  a  typical  fore¬ 
cast  entry  for  a  sample  predictand,  temperature.  Cate¬ 
gories  associated  with  a  temperature  forecast  are  listed, 
followed  by  the  forecast  entry: 

0  0  75  25  0. 

This  set  assigns  a  probability  of  occurrence  of  0%  to 
categories  1  and  2,  [-51  to  -19°C)  and  [-19  to  ~3°C); 

75%°  probability  of  occurrence  to  category  3  [~3  to  3°c)» 

25%  to  category  4  [3  to  11°C);  and  0%  to  category  5 
[ll  to  56°C]. 

The  weather  type  forecast  would  be  entered  in  this 

form : 

47  90 

which  indicates  weather  type  47,  light  freezing  drizzle, 
has  a  forecasted  probability  of  occurrence  of  90  percent. 

The  third  and  final  part  of  the  forecast  entry  dialog 
consists  of  the  questions  and  answers  which  build  a  complete 
forecast.  Fig.  D-tla  shows  a  sample  forecast  question. 

The  display  lists  forecast  station,  length  of  time  for 
which  the  forecast  is  being  made,  forecast  time  (initial 
time  plus  forecast  length),  and  predictand  being  forecast. 
This  information  is  followed  by  a  listing  of  the  categories 
for  that  predictand.  The  forecaster  then  enters  a  set  of 
probabilities  which  describes  his  forecast.  The  forecast 
entry  routine  checks  this  entry  to  assure  that  it  meets 
the  restrictions.  If  it  doer  not,  an  error  mess, age  is  dis- 
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played  which  identifies  the  restriction  violated  and  the 
question  is  redisplayed.  If  it  meets  the  restrictions 
the  data  are  entered  into  the  forecast  buffer  and  the  next 
question  is  asked.  This  process  continues  until  an  entire 
forecast  has  been  entered.  An  entire  forecast  consists  of 
nine  sets  of  questions  and  entries,  a  forecast  entry  for 
predictions  of  1,  3,  and  6  hours  for  each  of  the  three  pre- 
dictands.  When  the  entire  forecast  has  been  entered  it  is 
written  to  a  disk  file  known  as  the  forecast  file  and  the 
forecast  entry  dialog  is  terminated. 

C.  INTERROGATION  MODULE 

The  function  of  the  interrogation  module  is  to 
allow  the  forecaster  to  review  his  past  forecasts  and  his 
forecast  skill  scores. 

In  reviewing  past  forecasts  the  forecaster  is  given 
the  option  of  listing  any  or  all  of  his  five  most  recent 
forecasts.  Because  of  space  limitations  on  the  forecast 
file,  these  are  the  only  ones  saved.  After  he  enters  his 
selection,  a  hard  copy  of  the  forecast(s)  to  be  reviewed  is 
produced  on  the  line  printer  and  the  interrogation  dialog 
is  terminated. 

A  sample  output  is  shown  in  Fig.  D-12.  The  output 
header  lists  the  forecaster’s  name,  three  letter  ID  of  the 
forecast  station,  and  the  time  the  forecast  was  made  in 
the  form  MMDDHH  where  MM  is  the  month,  DD  is  the  day,  and 
HH  is  the  hour.  The  remainder  of  the  listing  is  divided 
into  nine  sections,  each  consisting  of  the  pertinent  data 
for  one  forecast  length  in  one  predictand.  Predict and  and 
forecast  length  begin  the  section,  followed  by  the  value  of 
the  predictand  at  the  time  the  forecast  was  made  (initial 
value).  The  value  of  the  predictand  at  the  verification 
time  (verification  value)  is  then  listed  if  it  is  avail¬ 
able.  The  probabilities  assigned  to  each  category  by  the 
forecaster  complete  the  data  in  one  section.  The  complete 
listing  provides  the  forecaster  with  information  to  determine 
in  a  subjective  manner  whether  he  made  a  correct  or  incor- 
rert  forecast. 
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***  AFGL  MCIOAS  ***  PROGRAM; 

REVIEW  FORECASTS  :  FORECASTER  MIKE  E  NIEOZIELSKI 


STATION:  BL V 

INITIAL  TIME:  13020 

predictand;  t 

FORECAST  LENGTH: 

1  HOURS 

INITIAL  VALUE  : 

-7 

OEGREES  C  VERIFICATION  VALUE 

CATEGORY  «  : 

1 

2  3  4 

3 

PROBABILITY  : 

0 

100  0  0 

0 

predictand:  T 

FORECAST  LENGTH: 

3  HOURS 

INITIAL  VALUE: 

-7 

DEGREES  C  VERIFICATION 

value: 

CATEGORY  *  : 

1 

2  3  4 

5 

PROBABILITY  * 

0 

100  0  0 

0 

PREDICTAND:  T 

FORECAST  LENGTH: 

6  HOURS 

INITIAL  VALUE: 

-7 

OEGREES  C  VERIFICATION 

value: 

CATEGORY  #  *• 

1 

2  3  4 

5 

PROBABILITY  : 

0 

100  0  0 

0 

PREDICTAND:  VIS 

FORECAST  LENGTH: 

1  HOURS 

INITIAL  VALUE  : 

1 

.20  MILES  VERIFICATION 

value 

CATEGORY  n  : 

1 

2  3  4 

5 

PROBABILITY  : 

0 

0  0  100 

0 

predictand:  VIS 

FORECAST  LENGTH: 

3  HOURS 

INITIAL  VALUE.* 

1. 

20  MILES  VERIFICATION  ' 

value : 

CATEGORY  #  : 

1 

2  3  4 

3 

PROBABILITY  : 

0 

0  0  O0 

20 

predictand:  vis 

FORECAST  LENGTH: 

6  HOURS 

INITIAL  VALUE: 

1 . 

20  MILES  VtR  I F I  CAT  I  ON  ’ 

value: 

CATEGORY  W  S 

1 

2  3  4 

5 

probability  : 

0 

0  0  20 

60 

PREOICTANO:  WX 

FORECAST  LENGTH: 

1  HOURS 

INITIAL  VALUE: 

s- 

VERIFICATION 

value:  non 

FORECASTED  VALUE 

• 

• 

S- 

CATEGORY  W  ' 

1 

2  3  4 

5 

PROBABILITY  :  100 

0 

PREDICTAND:  WX 

FORECAST  LENGTH: 

3  HOURS 

INITIAL  VALUE: 

s- 

VERIFICATION 

value:  S- 

FORECASTED  VALUE 

• 

• 

S- 

CATEGORY  W  : 

1 

2  3  4 

6 

PROBABILITY  ’•  100 

0 

PREDICTAND:  WX 

FORECAST  LENGTH: 

6  HOURS 

INITIAL  VALUE: 

s- 

VERIFICATION 

value:  s- 

FORECASTED  VALUE 

• 

• 

S- 

CATEGORY  *  : 

1 

2  3  4 

5 

PROBABILITY  :  100 

0 

ENO  OF  REVIEW  OF  FORECASTS 


FORCST 


S  MISSING 


-7  DEGREES  C 


-8  DEGREES  C 


MISSING 


.60  miles 


2.20  MILES 


Figure  D-12.  Sample  output  of  review  of  past  forecasts. 


The  review  skill  scores  option  allows  a  more  objective 
method  for  rating  forecasting  ability.  The  forecaster  se¬ 
lects  the  number  of  predictands  he  wants  to  review  and  then 
enters  the  predictands.  A  hard  copy  of  his  skill  scores 
for  those  predictands  is  then  produced  on  the  line  printer 
and  the  interrogation  dialog  is  terminated. 

A  sample  output  is  shown  in  Fig.  D-13-  The  output  is 
divided  into  sections  according  to  predictands  and  forecast 
lengths  similar  to  the  review  of  past  forecast  output. 

Each  section  is  identified  by  predictand  and  forecast  length. 
This  is  followed  by  the  number  of  forecasts  made  for  that 
predictand  and  forecast  length,  and  the  actual  skill  scores. 
These  scores  are  a  P-score,  cumulative  skill  score,  Allen 
utility  score,  percent  correct  skill  score,  and  Air  Weather 
Service  skill  score.  Biases  applied  to  each  category  when 
calculating  these  scores  are  then  listed,  completing  one 
section  of  data.  The  forecaster  can  use  these  scores  to 
chart  improvement  or  decline  in  his  forecasting  ability. 

D.  VERIFICATION  MODULE 

The  fourth  component  of  MFF  is  the  verification 
module,  whose  function  is  to  verify  the  forecasts  on  the 
forecast  file  and  calculate  individual  forecaster  skill 
scores.  This  module  is  invisible  to  the  user  since  once 
initiated  it  requires  no  further  terminal  inputs  to  perform 
its  function. 

When  initiated  the  verification  module  searches  the 
forecast  file  for  forecasts  which  have  not  been  verified. 
Encountering  such  a  forecast,  the  routine  retrieves  the 
actual  values  of  the  predictands  at  the  forecast  times  and 
incorporates  those  data  into  the  sector  of  the  forecast  file 
containing  the  unverified  forecast.  This  process  is  con¬ 
tinued  until  all  forecasts  on  the  file  have  been  verified. 
Then  the  probabilities  which  comprise  each  individual's 
forecast  are  compared  with  the  verification  value  or  a 
verified  persistence  forecast,  to  calculate  a  P-score,  Air 
Weather  Service  skill  score,  Allen  utility  score,  and  a 
percent  correct  skill  score.  A  cumulative  skill  score  for 


»»»  AFGL  MCIDAS  **«  PROGRAM.* 

SKILL  SCORE  REVIEW  :  FORECASTER  MIKE  E  NIEOZIELSKI 


predictand:  WX  forecast  length:  1  HOURS 

NUMRER  OF  FORECASTS  MADE:  2 

p-scope:  0.00000 

CUM  SKILL  SCORE:  0.00000 

ALLEN  UTILITY  SCORE:  1.000 

PER-CENT  CORRECT:  lOO.OOOOO 
AWS  SKILL  SCORE  UNDEFINED 
BIAS:  CATEGORY  1  :  -.02500 

BIAS:  CATEGORY  2  :  0.00000 


PREDICT  AND :  WX  FORECAST  LENGTH*.  3  HOURS 

NUMBER  OF  FORECASTS  MADE*.  2 

P- SCOPE:  1.00000 

CUM  SKILL  SCORE:  .50000 

ALLEN  UTILITY  SCORE:  .500 

PER-CENT  CORRECT:  50.00000 
AWS  SKILL  SCORE  UNDEFINED 
BIAS.*  CATEGORY  1  :  .95000 

BIAS:  CATEGORY  2  :  -.95000 


PREDICTAND:  WX  FORECAST  LENGTH:  6  HOURS 

NUMBER  OF  FORECASTS  MADE:  2 

P-SCOPE:  0.00000 

CUM  SKILL  SCORE:  0.00000 

ALLEN  UTILITY  SCORE:  1.000 

PER-CENT  CORRECT.*  1  00.00000 
AWS  SKILL  SCORE:  1.00000 
BIAS:  CATEGORY  1  :  -.02500 

BIAS:  CATEGORY  2  :  0.00000 


FORCST 


Figure  D-13- 


Sample  output  of  review  of  skill  scores 


all  the  individual's  forecasts  for  each  predictand  and  fore 
cast  length  is  also  calculated,  and  the  category  biases  for 
doing  these  calculations  are  updated.  This  process  is  re¬ 
peated  for  each  forecaster  ID  on  the  forecast  file.  Then  a 
message  indicating  completion  of  verification  is  displayed 
and  control  is  returned  to  the  beginning  of  the  MFF  program. 
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E.  IMPROVEMENT  OF  BRIGHTNESS  ANALYSIS  TECHNIQUE 


I.  INTRODUCTION 

Use  of  geostationary  satellite  data  to  provide  infor 
mation  regarding  cloud  cover  and  significant  weather  re¬ 
quires  the  development  of  algorithms  to  discriminate  be¬ 
tween  the  earth's  surface  and  cloud  type  and  amount  in 
the  automatic  processing  of  satellite  imagery.  SASC 
advances  in  this  area  are  the  subject  of  this  section. 

II.  DATA  SAMPLE 

Data  to  be  analyzed  by  the  processing  software  con¬ 
sist  of  hourly  visible  and  infrared  imagery  for  the 
northeastern  U.  S.  as  received  from  two  geostationary 
satellites  (GOES  I  and  GOES  II).  The  region  covered 
corresponds  to  5°°  visible  scan  lines,  each  consisting 
of  765  picture  elements  or  pixels  of  one-mile  resolution 
brightness  counts.  The  IR  imagery  consists  of  125  scan 
lines  each  containing  896  pixels  of  four-mile  resolution 
thermal  counts.  The  IR  sensor  samples  every  two  miles 
(at  the  subpoirt)  along  the  scan  line,  resulting  in  twice 
as  many  pixels  in  the  east-west  direction  as  in  the  north 
south  direction. 

One  hour's  data  are  preceded  by  a  header  record  con¬ 
taining  Greenwich  Mean  Time,  first  and  last  line  numbers 
of  the  IR  rows  recorded,  and  the  number  of  the  first  ele¬ 
ment  in  the  visible  lines  to  be  recorded.  These  row  and 
element  numbers  are  based  on  their  location  in  a  full 
disk  picture  of  the  earth.  The  data  sample  is  divided 
into  1? 5  sets  consisting  of  four  visible  data  lines,  one 
per  record,  and  one  record  of  IR  data.  Each  visible 
scan  line  record  is  preceded  by  a  standard  identification 
word.  The  IR  record  also  contains  128  words  of  data 
which  describe  the  Greenwich  Mean  Time  at  which  the  scan 
occurred,  Julian  date,  line  number,  and  other  information 
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not  relevant  here.  Each  count  value  is  an  8-bit  byte 
and  is  recorded  as  part  of  a  24-bit  word.  Each  tape  con¬ 
tains  a  maximum  of  24  images,  usually  for  the  six  to 
eight  hours  when  the  sun  provides  sufficient  illumination; 
i.e.,  local  afternoon. 

Location  of  any  point  in  the  longitude-latitude  co¬ 
ordinate  system  can  be  expressed  in  the  line-element 
coordinate  system  through  the  use  of  a  navigation  trans¬ 
form  based  on  the  full  earth  disk  image  scanned  by  the 
satellite  and  known  orbital  parameters.  Line  and  ele¬ 
ment  numbers  of  four  standard  locations  within  the  data 
area  (46°N,  83°W,  and  the  FAA  reporting  stations  at  Bed¬ 
ford,  Mass.,  Block  Island,  R.I.,  and  Dulles  International 
Airport,  Va.)  are  used  as  input  data  for  reference 
points  to  locate  any  other  point  in  the  data  area. 

III.  SOFTWARE  DESCRIPTION  AND  DEVELOPMENT 

Data  processing  software  consists  of  four  parts: 

1)  station  locating  routine  (program  CENTER)  which  re¬ 
trieves  imagery  data  for  a  particular  area  from  the  data 
tape;  2)  bit  manipulating  routine  (subroutine  EIBYT)  to 
convert  the  24-bit  word  format  of  the  data  tape  to  the 
60-bit  word  format  of  the  analyzing  hardware;  3)  data 
normalizing  routine  (subroutine  NORMAL)  to  correct  mis¬ 
sing  or  garbled  data;  and  4)  data  analyzing  routine  (sub¬ 
routine  BRIGHT). 

CENTER  calculates  the  row  and  element  number  of  the 
station  for  each  hour  to  be  analyzed.  The  data  buffer  is 
then  filled  with  the  brightness  value  of  each  pixel  in 
an  array  of  7,200  elements  (120  rows  by  60  elements) 
using  the  station  row  and  element  as  the  center  and  the 
line  number  associated  with  each  row  of  data.  The  large 
size  of  this  array  insures  an  adequate  sample  for  the  nor¬ 
malizing  routine  to  function  properly.  While  each  row  of 
data  is  being  entered  into  the  array,  EIBYT  is  converting 
the  24-bit  data  words  containing  three  8-bit  brightness 
values  into  60-bit  words  containing  one  brightness  value. 
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This  process  is  necessary  to  resolve  the  data  structure 
differences  between  the  archiving  hardware  (McIDAS)  and 
the  analyzing  hardware  (CDC  6600).  The  data  buffer  is 
then  passed  to  NORMAL  and  BRIGHT  for  normalization  and 
analysis.  Upon  completion  of  analysis  the  data  tape  is 
advanced  to  the  next  hour's  data,  the  reference  point  re¬ 
located,  and  the  process  is  repeated,  until  all  hours  on 
the  tape  have  been  analyzed. 

In  an  effort  to  increase  the  amount  of  analysis 
done  during  each  run  of  the  data  tape,  CENTER  was  modified 
to  calculate  the  line  and  element  number  of  up  to  five 
stations.  The  number  of  lines  between  the  station  being 
analyzed  and  the  next  station  on  the  list  can  be  deter¬ 
mined  algebraically.  This  allows  the  tape  to  be  advanced 
or  rewound  the  appropriate  number  of  lines  so  as  to  be 
positioned  at  the  beginning  of  the  120  lines  of  data  cen¬ 
tered  on  the  new  station.  In  this  fashion  several  stations 
can  be  analyzed  before  the  tape  is  advanced  to  the  next 
hour,  thus  increasing  the  amount  of  analyzed  data  per  run 
as  much  as  five  times . 

Since  some  data  are  inevitably  lost  or  garbled  in 
transmission  from  the  satellite,  NORMAL  was  developed  to 
normalize  or  correct  them.  Normalization  begins  by  cal¬ 
culating  the  average  of  each  of  the  120  lines  of  data. 

If  a  line  average  differs  by  more  than  20  from  the  pre¬ 
ceding  line  average,  that  line  is  flagged  to  indicate  that 
it  contains  garbled  or  noisy  data.  To  prevent  comparing 
a  line  average  to  the  average  of  a  garbled  line  of  data, 
the  average  of  a  noisy  line  is  replaced  by  the  average  of 
the  last  good  line  of  data  encountered.  After  all  lines 
in  the  sample  have  been  tested,  noisy  lines  are  replaced 
element  for  element  by  the  average  of  corresponding  ele¬ 
ments  in  the  two  adjacent  good  lines  rf  data. 

The  final  step  involves  isolating  the  27  lines  of 
data  which  will  be  processed  by  the  analyzing  routine. 

This  is  done  by  calculating  the  line  number  of  the  thir¬ 
teenth  line  of  data  before  the  station  row.  This  line  and 
the  next  26  lines  of  data  are  then  copied  onto  a  temporary 
disk  file  to  be  passed  to  BRIGHT  for  analysis. 
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Several  modifications  were  made  to  this  routine  to 
reduce  processing  time  and  increase  data  integrity.  First, 
the  line  numbers  of  the  27  lines  of  data  to  be  analyzed 
are  checked  to  insure  that  they  are  consecutive  lines. 

The  original  normalizing  process  ignored  missing  lines  by 
merely  taking  the  first  2?  lines  of  date  it  encountered 
which  had  line  numbers  greater  than  or  equal  to  the  initial 
row  number  of  the  sample  set  to  be  analyzed.  The  modified 
routine  can  detect  non-consecutive  line  numbers  and  calcu¬ 
late  the  number  of  imbedded  blank  lines.  These  lines  are 
filled  with  dummy  brightness  values  and  flagged  as  being 
noisy.  The  process  of  isolating  the  2?  line  set  was 
shifted  to  occur  before  the  replacing  of  noisy  lines 
takes  place.  This  allows  blank  lines  to  be  filled  with 
normalized  data  and  also  allows  the  routine  to  replace 
noisy  lines  within  the  data  sample  to  be  analyzed  without 
having  to  do  it  for  the  entire  120  line  normalizing  sample. 
To  further  reduce  normalizing  time,  the  use  of  a  scratch 
disk  file  to  hold  the  normalized  data  and  its  attendant 
input-output  was  eliminated  by  placing  only  the  27  ele¬ 
ments  of  each  of  the  27  lines  of  data  to  be  analyzed  into 
a  common  block  of  storage  within  the  program.  This  pro¬ 
vides  instant  program  access  without  the  need  for  time- 
consuming  reading  from  or  writing  onto  a  scratch  file. 

After  the  data  have  been  normalized  and  isolated 
the  analysis  routine  BRIGHT  begins  processing.  The  27  x 
27  array  is  divided  into  nine  9x9  blocks,  as  shown  in 
Fig.  E-la,  and  each  of  these  blocks  is  further  subdivided 
into  3x3,  5x5,  7  x  7,  and  9x9  pixel  boxes  (see  Fig. 
E-lb)  for  analysis  purposes.  Average  brightness,  stan¬ 
dard  deviation,  and  the  range  of  values  for  each  box  are 
calculated  and  the  maximum  and  minimum  brightness  value 
in  each  are  determined.  The  brightness  value  of  the  cen¬ 
tral  pixel  of  these  nested  boxes  is  noted  and  the  gradient 
from  each  edge  to  the  center  is  calculated  for  all  four 
boxes . 

These  analyses  are  performed  for  each  of  the  nine 
blocks  and  for  the  entire  27  x  27  array.  A  sample  output 
of  tee  analysis  is  shown  in  Fig.  K-2.  The  first  two 
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Figure  E-la.  Divisions  of  the  2?x27  array.  Dimen¬ 
sions  are  in  pixels. 


Figure  E-lb.  Nested  boxes  within  a  9*9  block. 

Shaded  area  is  central  pixel. 
Dimensions  are  in  pixels. 
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Figure  E-2.  Sample  brightness  analysis  output.  Values  are  visual 
count  values  based  on  a  256-count  range. 


lines  of  output  give  the  Julian  day  and  Greenwich  Mean 
Time  of  the  data  being  analyzed,  station  name,  its  row  and 
element  numbers,  and  a  fine-tuned  navigation  correction 
factor  which  aligns  the  image  coordinates  with  the  earth 
surface  coordinates.  The  label  which  precedes  the  hourly 
data  on  the  tape  is  also  listed.  Analysis  of  the  nine 
blocks  is  arranged  as  shown,  followed  by  analysis  of  the 
27  x  27  array.  A  final  line  indicates  the  number  of 
blank  lines  of  data  in  the  analysis  sample.  This  is  in¬ 
cluded  to  provide  a  basis  for  accepting  or  rejecting  the 
analysis.  All  of  these  output  values  are  also  sorted  by 
station  and  written  to  separate  permanent  files  for  future 
reference  and  further  analysis. 

To  augment  this  brightness  value  analysis  and  elimin¬ 
ate  data  dependency  on  varying  solar  elevation  angles, 
routine  ALBEDO  was  written  to  convert  brightness  analy¬ 
sis  results  to  their  corresponding  albedo  values.  It 
was  decided  to  perform  this  task  using  a  separate  routine 
to  keep  CENTER  running  time  to  an  absolute  minimum.  The 
maximum,  minimum,  central,  and  range  of  values  for  bright¬ 
ness  data  can  be  converted  to  albedo  values  by  using  the 
formula : 

A  =  (B/C)2  (1/cos  Z) 

where  A  is  albedo  value,  B  is  corresponding  brightness 
value,  C  is  a  constant  equal  to  for  GOES  I  and  255 
for  GOES  II  data,  and  Z  is  solar  zenith  angle.  Z  is  cal¬ 
culated  by  a  separate  routine  using  Julian  day,  Green¬ 
wich  Mean  Time,  and  longitude  and  latitude  of  the  station. 

The  average  and  standard  deviation  conversion  was 
more  involved  since.’ 

A  =]>](B/C)2  (|j)(  1/cos  Z) 


and 


A 


(B/C)  ( 1/ cos  Z) 


L/N  -  A£ 


and  B  and  B  are  not  directly  calculated  by  BRIGHT 

t—J  V— >  2 

and  included  in  the  output  data.  >  B  '  can  however  be 
determined  from  the  equation: 

0  B  =  [Eb2/N  "  ^ 
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since  and  B  are  known  values.  N  can  be  calculated 

n 

from  the  analysis  box  size.  Only  a ^  cannot  be  calculated 
from  data  output  by  BRIGHT.  This  was  corrected  by  add i nr 
B  to  the  values  calculated  in  that  routine  and  saved 
on  the  output  file.  Gradient  analyses  were  omitted  from 
the  albedo  analysis  because  no  new  information  could  be 
obtained  from  them  in  albedo  form. 

Output  format  for  albedo  analysis  ( Fig .  E-3)  is 
similar  to  the  brightness  analysis.  Use  of  brightness 
and  albedo  data  provides  a  basis  for  discriminating  be¬ 
tween  surface  brightness  and  clouds. 

Efforts  to  include  infrared  data  analysis  in  CENTER 
have  begun.  Differences  between  one-mile  resolution 
brightness  data  and  four-mile  IR  data  have  been  recon¬ 
ciled  by  dividing  each  IR  pixel  into  8  pixels  of  the 
same  value.  These  expanded  IR  data  are  then  aligned 
with  the  brightness  data  to  produce  a  one-to-one  corres¬ 
pondence  betweeen  brightness  or  visible  data  pixels  and 
IR  pixels.  Brightness  analysis  software  can  then  be  used 
on  IR  data.  Difficulties  in  normalizing  IR  data  have 
thus  far  prevented  sustained  analysis. 

IV.  CONCLUSIONS 

Improvements  in  program  CENTER  have  increased 
efficiency  of  execution  time  and  reliability  of  results. 
Run  time  per  station  analyzed  has  been  reduced  62  per¬ 
cent  and  field  length  requirement  has  been  halved,  re¬ 
sulting  in  a  32  percent  cost  saving  per  station  analyzed. 
Incomplete  data  samples  can  be  identified  and  rejected 
on  the  basis  of  the  number  of  blank  lines  of  data,  re¬ 
moving  a  source  of  erroneous  results  in  later  analyses. 
The  addition  of  albedo  analysis  has  increased  the  data 
base  for  ground/cloud  discrimination. 
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